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Abstract. Primary inversions of accurately measured solar oscillation frequencies coupled with the equations of 
thermal equilibrium and other input physics, enable us to infer the temperature and hydrogen abundance profiles 
inside the Sun. These profiles also help in setting constraints on the input physics as well as on heavy element 
abundance in the solar core. Using different treatments of plasma screening for nuclear reaction rates, limits on 
the cross-section of proton-proton nuclear reaction as a function of heavy element abundance in the solar core are 
obtained and an upper limit on heavy element abundance in the solar core is also derived from these results. 
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1. Introduction 

The precisely measured frequencies of solar oscillations 
have been used to probe the solar interior. The primary 
inversions of these observed frequencies yield the sound 
speed and density profiles inside the Sun. In order to in- 
fer the temperature and chemical composition profiles, 
we also need to know the input physics such as opaci- 
ties, equation of state and n uclear energ y gene ration rates 
(Gough & Kosovichev 1988 ; Kosovic hev 1996 ; Shibahashi 
& Takata ^9961 Takata & Shibahashi Il998l Antia & Chitre 



1998| ). In all these works the heavy element abundance 
profile is assumed to be known; attempts to determine 
heavy element abundance profile from helioseismic data 
have not been particularl y succ essful (Antia & Chitre 
199£; Takata & Shibahashi 2001 ) as the resulting inverse 
problem becomes extremely ill-conditioned. Fukugita & 
Hata ( 199^ ) have obtained limits on heavy element abun- 
dance in the solar core using observed solar neutrino 
fluxes. It would be interesting to enquire if such limits 
can be independently obtained from helioseismic data. 

In general, the computed luminosity in a seismically 
computed solar model is not expected to match the ob- 
served solar luminosity. By applying the observed lu- 
minosity constraint it is possible to constrain the in- 
put physics, particularly, the cross-section of proton- 
proton (pp) nuclear reaction. Antia & Chitre ( |1998| ) es- 
timated this cross-section to be Sn = (4.15 ± 0.25) x 
10~^^ MeV barns. Similar values have been obtained by 



comparing the computed solar models with helioseismic 
data (DeglTnnocenti, Fiorentini & Ricci 1998; Schlattl, 



Bonanno & Paterno 1999). The main source of error in 



these estimates is the uncertainty in the Z profile and, 
therefore, Antia & Chitre ( 1999| ) attempted to find the 
pp reaction rate as a function of Z in the solar core. In 
all these works the plasma screening of nuclear reaction 
cross-sections was calculated using intermediate screen- 



ing formulation of Graboske et al. (1973). The treatment 
of screening in stellar nuclear reaction rates is not yet 
adequately understood (Dzitko et al. 1995; Gruzinov & 
Bahcall|l99|). Wilets et al. ( pOOOj ) have done a sophisti- 
cated treatment of plasma screening and compared their 
results with earlier prescriptions. Their results indicate 
that for the solar core the intermediate screening treat- 
ment due to Mit ler (1977 ) is better than that due to 
Graboske et al. ( 1973 ). It would thus be interesting to 



study the effect of different treatment of plasma screen- 
ing on the helioseismically estimated pp reaction cross- 
section. 



Antia & Chitre (1999) included the effect of heavy ele- 
ment abundance Z only on the opacity of the solar mate- 
rial. If we make the reasonable assumption that the abun- 
dances of C, N, O also increase with Z, then the CNO cycle 
will become more effective in contributing to the nuclear 
energy generation in the solar core. At normally accepted 
values of Z it is estimated that less than 2% of energy 
generated in the central region is produced by the CNO 
cycle (Bahcall et al. 2001). But if Z value is increased, 
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this proportion will clearly increase and consequently, the 
pp reaction rate needs to be reduced to maintain the ob- 
served solar luminosity. In this work we demonstrate that 
this effect can be exploited to set an upper limit on Z in 
the solar core. 
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2. The technique 

With the use of accurately measured p-mode frequen- 
cies from the first year of operation of MDI instrument 



4.5 



(Rhodes et al. 1997) we infer the sound speed and density 
profiles adopting a Regularised Least Squares technique 
(Antia |1996 ). The determination of thermal and chemi- 
cal composition profiles, however, necessitates the use of 
equations of thermal equilibrium (Antia & Chitre 1998) 
with the supplementary input of the equation of state, 
opacities and nuclear energy generation rates, provided we 
have a knowledge of the heavy element abundance pro- 
file inside the Sun. There is no guarantee for the resul- 
tant seismic model to yield the observed solar luminosity, 
Lq = 3.846 X 10'^'^ ergs s~^, unless we adjust the nuclear 
reaction rates slightly. This turns out to provide a valuable 
tool for constraining the nuclear reaction rates, in partic- 
ular, the cross-section of the pp nuclear reaction which 
has not been measured in the laboratory because of the 
slow weak interaction rate. The pp reaction cross-section 
has only been calculated theoretically and it would be in- 
structive to test the validity of calculated results using the 
helioseismic constraints. 

There is an uncertainty of about 2% in evaluating the 
luminosity of seismic models because of possible errors 
in primary inversions (including an estimate of system- 
atic errors), solar radius, equation of state, opacity, nu- 
clear reaction rates for other reactions. Possible error due 
to uncertainties in treatment of plasma screening of nu- 
clear reactions is not included in this error estimate. This 
can be estimated by using different treatments of plasma 
screening. The uncertainty arising from errors in Z pro - 
files is, however, much larger (cf., Antia & Chitre 1998| ). 
We, therefore, estimate for each Z profile, the range of 
cross-section of pp nuclear reaction, which reproduces the 
luminosity to within 2% of the observed value. The com- 
puted luminosity in the seismic models will, of course, de- 
pend on the metallicity, Zc in the solar core; we therefore, 
attempt to delineate the region in the Sn-Zc plane which 
gives the correct solar luminosity. In order to obtain the 
thermal structure we adopt the OPAL opacities (Iglesias 
& Rogers 1996), the OPAL equation of state (Rogers, 
Swenson & Iglesias 199t) and the nuclear reaction rates 
from Adelberger et al. (1998). The plasma screening ef- 
fects are calculated using either weak screening (Salpeter 



1954) or intermediate screening (Mitler 1977). The effect 
of variation of Zc on nuclear energy generation rate is 
incorporated in the computations by assuming the abun- 
dances of all heavy elements to increase in the same ratio 
for obtaining the abundances of C, N, O. 

3. Results 

With the help of the inverted profiles for sound speed and 
density, along with a homogeneous Z profile, covering a 
wide range of Z values, we first calculate a seismic model 
by employing the equations of thermal equilibrium. For 
each central value of Z we estimate the range of cross- 
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Fig. 1. The region in ZcSu plane that is consistent with 
helioseismic data is marked by horizontal shading. The 
continuous line defines the values obtained using interme- 
diate plasma screening, where the seismic model matches 
the observed solar luminosity. The point with error bars 
shows the current best estimates for Zc and The 
short-dashed line shows the central value obtained using 
weak screening, while the long-dashed line shows the cen- 
tral value obtained using intermediate screening and also 
including the effect of Z on nuclear energy generation. 
The two horizontal lines mark the limits of theoretically 
calculated values of ^n. 



section of pp nuclear reaction, which reproduces the lumi- 
nosity to within 2% of the observed value. The numerical 
results are shown in Fig. 1, which delineates the region in 
Zc^Sii plane that is consistent with helioseismic and lumi- 
nosity constraints. This figure shows the results obtained 
using intermediate screening (Mitler 1977). For compar- 
ison the central value obtained using weak screening is 
also shown in this figure. In order to compare these re- 
sults with the earlier work of Antia & Chitre (1999), in 
these calculations the effect of Z on nuclear energy gen- 
eration rate is not included. Thus the difference is mainly 
due to treatment of plasma screening. It is clear that with 
these treatments of screening the curren t estim ate of pp 
reaction cross-section (Adelberger et al. 1998| ) is consis- 
tent with helioseismic estimate within the expected error 
bars, though the best value is about 1.5% higher than 
the theoretical estimates. Further, there is not much dif- 
ference between results obtained using the two different 
treatments of plasma screening. In general, the estimated 
value of Sii using weak screening is smaller than that ob- 
tained using intermediate screening by about 1%. 

With a view to study the effect of Z on energy gener- 
ation through CNO cycle, we repeat the calculations by 
varying the abundances of all heavy elements in the same 
proportion as Z in the nuclear reaction network. These 
results are shown by long-dashed line in Fig. 1. We do not 
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Fig. 2. The seismically estimated value of Sn as a func- 
tion of Zc the heavy element abundance at the centre for 
different Z profiles. The solid line shows the value obtained 
with homogeneous Z profile and is the same as the long- 
dashed line in Fig. |^. The short-dashed and long-dashed 
lines show the results obtained using Z profile given by 
Eq. |l| and respectively. The dotted line shows the re- 
sults using Z profile of Eq. |I|, but with the reaction rate 
of + p reaction reduced by 40%. The point with error 
bars shows the current best estimates for Zc and Sn. 



expect much difference at low values of Zc, since CNO 
cycle produces a tiny fraction of energy in the Sun. But 
with augmented Zc, this fraction becomes larger due to 
increased Z as well as due to the resulting increase in the 
temperature, since the CNO reactions are more sensitive 
to temperature. As a consequence, in order to maintain the 
solar luminosity constraint we need to reduce ^n, with the 
result that around Zc = 0.035 the estimated value of Su 
decreases very sharply. This is because a significant frac- 
tion of the luminosity is then accounted for by the CNO 
reactions with seismically inferred temperature and com- 
position profiles. For Zc = 0.035, CNO cycle accounts for 
about 8% of solar luminosity and this fraction, in fact, in- 
creases rapidly with Zc- Thus, we can consider this as the 
upper limit on Zc, which is somewhat larger than that 
inferred by Fukugita and Hata ( f 998] ). A lower limit on 
Zc can also be obtained from this figure if we accept the 
range of theoretically computed values of 5ii (Bahcall & 
Pinsonneault 1995; Turck-Chieze & Lopes 1993| ) as repre- 
senting the acceptable range of ^n. Assuming an upper 
limit of 4.3 x 10~^^ MeV barns for 5*11 we get a lower limit 
of Zc = 0.006, which is again comparable to the lower limit 
obtained by Fukugita and Hata ( 1998| ). If we adopt a lower 
limit of 3.8 x 10^^^ MeV barns for 6*11, the resulting up- 
per limit on Zc = 0.032. The exact limits on Zc will of 
course, depend on the assumed limits on Su, but because 
of the sharp fall in the long-dashed curve in Fig. 1 , the up- 
per limit is not particularly sensitive to the assumed lower 



Table 1. Neutrino fluxes in a seismic model using the Z 
profile from model NO of Brun et al. (|2002). 



Source Flux (cm s ) 

(6.05 ±0.06) X 10'" 
(1.42 ±0.02) X 10* 



PP 
pep 
hep 
^Be 
«B 

13n 

Total CI 

Total Ga 



2.09 X 10^ 
(4.76 ± 0.48) X 10^ 
(4.83 ± 0.88) X 10*^ 
(5.27 ±0.79) X 10* 
(4.48 ± 0.80) X 10* 
(2.76 ± 0.30) X 10*^ 
7.27 ± 1.14 SNU 
127.8 ± 7.2 SNU 



limit on 5*11 and any value Zc > 0.035 would certainly be 
difficult to reconcile with seismic models. 

In the foregoing calculations we have adopted a ho- 
mogeneous Z profile, which is perhaps not too realistic. 
Since most of the energy generation occurs in the solar 
core (r < O.25i?0) our results are unaffected by the Z 
profile in the outer radiative region. To demonstrate this 
we also use the following Z profiles : 



Z = 



and 



0.463 



{r/RQ -0.25) 



if r < 0.25i?o, 
otherwise, 



Z = Zc + iZ,- Zc)r/{0.713RQ), 



(1) 



(2) 



where Zg = 0.018 is the known value of Z at the solar 
surface. The seismically estimated value of 5*11 for these 
Z profiles is compared with those for the homogeneous 
profile in Fig. ||. As to be expected, for Z profile given by 
Eq. |l| the results arc almost identical to those obtained 
with the homogeneous profile, while for the linear Z pro- 
file given by Eq. ^ the average value of Z in the solar core 
would be somewhat lower and as a result the estimated 
Su values are somewhat larger, but the difference is com- 
parable to the estimated errors. If instead of central value 
of Z we had used the value of Z at around O.li?© for 
the linear profile, the two curves would have been almost 
identical. Thus it is clear that exact form of the Z profile 
is not very important and the upper limit on Zc is only 
weakly dependent on the Z profile. 

Apart from uncertainties in Z profile the nuclear re- 
action rate for the CNO cycle reactions are also rather 
uncertain and these may affect the inferred upper limit. 
To estimate this uncertainty we repeat the calculations us- 
ing the profile given by Eq. ^ for Z, but with the reaction 
rate of the -^^ISl + p reaction reduced by 40%, which is the 
estimated uncertainty in this reaction rate (Adelberger et 
al. 1998D_ . These results are also shown by the dotted line 
in Fig. g. From the figure it is clear the uncertainties in 
CNO reaction rates will not particularly affect the upper 
limit on Zc in a significant way. 

Table 1 lists the neutrino fluxes in a seismic model ob- 
tained using the Z-profile of a model including tachocline 
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mixing of Brun et al. (2002), with the use of intermediate 
screening, for calculating the nuclear energy generation 
rates. In this case we need to increase 5*11 by 1.6% over the 
currently accepted value. As noted by Brun et al. ( 2002D , 
an increase in Sn leads to a better agreement between 
solar models and the seismically inferred sound speed and 
density profiles. Thus a small increase in Sn may still be 
required, but it is within the errors in helioseismic esti- 
mates and those in theoretically computed values. These 
neutrino fluxes are generally larger than those in seismic 
model of Antia & Chitre (1998), presumably because of re- 
duction in 5*1 1. The ^B neutrino flux in the seismic model 
is consistent with the current estimate of (5.09±0.65) x 10^ 
cm~^ s~^ using the neutral current channe l of S udbury 
Neutrino Observatory (SNO) (Ahmad et al. ^002|) . 



4. Discussion and Conclusions 

With the help of inverted sound speed and density pro- 
files, it is possible to infer the T, X profiles in the solar 
interior, provided the Z profile and the input physics are 
known. The resulting seismic models have the correct solar 
luminosity, provided the heavy element abundance Zc in 
the solar core and the cross-section for pp nuclear reaction 
are within the shaded region shown in Fig. 1. It appears 
that the currently accepted values of Zc or 6*11 need to 
be increased marginally to make them consistent with he- 
lioseismic constraints. The required increase is within the 
error estimates. The higher estimates for Sn obtained ear- 
lier were due to differences in treatment of plasma screen- 



ing. With the use of weak (Salpeter 1954) or intermediate 
screening due to Mitler (1977) the theoretically estimated 
value of 5*1 1 is in reasonable agreement with seismically es- 
timated value. W ith a Z profile in a standard solar model 
NO of Brun et al. (2002), the seismically estimated value of 
Sii is 4.07 X lO^'^^ MeV barns with the use of intermediate 
screening and 4.02 x 10~^^ MeV barns for weak screening 
used while calculating the nuclear energy generation rate. 

If the value of heavy element abundance in the solar 
core, Zc, is increased beyond 0.035 the CNO cycle gen- 
erates a good fraction of solar luminosity and it is not 
possible to get any consistent seismic model unless Sn 
is decreased substantially below the accepted value. This 
puts a clear upper limit on the heavy element abundance 
in the solar core, which is comparable to that indepen- 
dently obtained by Fukugita and Hata (1998). This upper 
limit is not very sensitive to the Z profile or to the uncer- 
tainties in the CNO reaction rates. Even if the additional 
heavy elements in solar core do not include CNO, the ef- 
fect of opacity alone will also put an upper limit on Zc, but 
in that case the limit will depend on the assumed lower 
limit on Sn. For currently accepted theoretical limits, the 
upper limit in this case turns out to be around 0.035. 

Note that the seismic model satisfies the normal stellar 
structure equations, though the inferred X profile may not 
match that given by an evolutionary solar model. Further, 
since the seismic model is confined to the radiative inte- 
rior, it will not be possible to match it to an acceptable 



convection zone model, unless the heavy element abun- 
dance at the top of the radiative region is close to the 
known surface value. Thus a relatively high value of Z in 
the solar core can be realised only if the Z profile has a 
significant gradient in the radiative region. From experi- 
ments with varying Z or Sn in evolutionary solar models 
also it is known that an increase in Z can be compensated 
for by a reduction in 5*1 1 to match the seismically inferred 
sound speed and density profiles (Brun et al. p002 ). Thus 
it is quite possible that a comparable upper limit on Zc 
may be obtained from these models if we restrict the range 
of S*!!, provided the initial Z profile at the zero age main 
sequence stage is not homogeneous. 
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